Dynamic Failure in Amorphous Solids via a Cavitation Instability 
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The understanding of dynamic failure in amorphous materials via the propagation of free bound- 
aries like cracks and voids must go beyond elasticity theory, since plasticity intervenes in a crucial 
and poorly understood manner near the moving free boundary. In this Letter we focus on failure 
via a cavitation instability in a radially-symmetric stressed material, set up the free boundary dy- 
namics taking both elasticity and visco-plasticity into account, using the recently proposed athermal 
Shear Transformation Zone theory. We demonstrate the existence (in amorphous systems) of fast 
cavitation modes accompanied by extensive plastic deformations and discuss the revealed physics. 



The principal difficulty in describing the mechanical 
failure of materials using elasticity theory is that this 
approach does not provide a law of motion for the free 
boundary that naturally occurs when a crack or a void 
is propagating. Recently, phenomenological phase-field 
models [l| were proposed in order to overcome this fun- 
damental problem. In reality the high stress concentra- 
tion near the moving boundary must be associated with 
some form of plastic deformations that are rather poorly 
understood in the context of amorphous materials; the 
consequences of these plastic deformations are usually 
referred to as the "process zone" whose properties one 
usually does not dare to probe too closely. The aim of 
this Letter is to do just the opposite 0. We set up the 
simplest (from the point of view of symmetries) example 
of free boundary dynamics, i.e. that of a circular cavity 
responding to radially symmetric stresses at infinity (or 
equivalently high pressure inside the cavity). This prob- 
lem has a rather long history 0, 0] ; the main contribution 
of our study is the elucidation of the role of amorphous 
plasticity and the detailed interaction between elasticity 
and plasticity, throughout the system and in particular 
near the free boundary where plastic dynamics is explic- 
itly described in terms of the athermal Shear Transforma- 
tion Zone theory (STZ) 0,0]. This theory automatically 
includes hardening and rate-dependent effects in addi- 
tion to a proper Eulerian description of the equations of 
motion which allows the discussion of inertial effects and 
of large deformations including accelerating and catas- 
trophic cavitation instability. 

Consider then an infinite 2D medium with a circular 
hole around the origin of radius R(t), loaded by a radially 
symmetric stress a°° at infinity. The symmetry dictates 
a radial velocity field v that is independent of the az- 
imuthal angle 8, i.e. v r (r,t) = v(r,t) and vg(r,t) = 0. 
This velocity field implies the components of the total 
rate of deformation tensor 
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In this Letter we restrict the elastic part of the deforma- 
tion to be small, allowing us to decompose the total rate 
of deformation tensor into its elastic and plastic parts. 



Denoting the material time derivative as 'D t =dt + vd r 



where the plastic part DfJ will be discussed below. The 
(linear) elastic part e|j in the present symmetry is deter- 
mined by the deviatoric stress s=see=—s rr =(crgg—cr rr )/2 
and the pressure p= — (agg + a rr )/2 according to 

4 = -p/2K-s/2p, 
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where <Xy is the stress tensor and K and p are the 2D bulk 
and shear moduli respectively. The plastic rate of de- 
formation is assumed traceless, corresponding to incom- 
pressible plasticity. In the present symmetry the plastic 
rate of deformation tensor has only one independent com- 
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and ([3]) in Eq. ([2]) one ends up with 
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At this point we introduce the equations of motion for 
the material density p and the velocity v in the form 

V tP = -pV ■ v , (6) 
P V t v = V-* = -d r (r 2 s) /r 2 - d rP , (7) 

where the boundary itself evolves according to R(t) = 
v(R,t). The boundary conditions are given by 



a rr (R,t) = -p(R,t) -s(R,t) = , 
oy r (oo,i) = — p(oo, t) — s(oo, i) = er° 



(8) 



As initial conditions for the dynamics we take the static 
linear elastic solution for a medium with the given stress 
at infinity and a hole of radius i?(0), assuming that the 
initial elastic response is much faster than the plastic 
deformation time scale, 



P(r,0) 



,(r,0) = a°°i? 2 (0)/r 2 , u(r,0) = 



(9) 
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Finally, we supply the constitutive relations determin- 
ing D pl . The STZ theory determines this object in 
terms of a low density A of shear transformation zones; 
these are localized regions embedded in the elastic matrix 
causing an irreversible deformation in response to shear 
stress. Under such a stress they are assumed to transform 
between two internal states, resulting in a plastic strain 
increment along one of the two principal axes in our sym- 
metry. Macroscopic plastic flow, denoted by D pl , results 
from the average flippings of STZ's. This flow, by itself, 
can also create and annihilate other STZ's. The bias of 
the populations of STZ between these two states, caused 
by the application of a deviatoric stress, is denoted by a 
tensor m which in our symmetry has only one indepen- 
dent component denoted by to. This field is a local order 
parameter acting as a "back stress" . The total density of 
STZ's is affected by a normalized effective temperature 
X, characterizing the state of configurational disorder of 
the material, via a Boltzmann-like factor. Details of the 
derivation of the equations can be found in Q . Here we 
just state the final equations where s had been normal- 
ized with the dynamic yield stress s y , which is a material 
parameter, 
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{\s\ - s a ) s c a exp(-C s a ) ds a . (14) 



The parameter r sets the basic time scale of plasticity. c 
and eo are parameters of the order of unity that in this 
work are taken as unity. The function C(s) quantifies 
the stress-dependence of the transition rate of individual 
STZ's between their internal states, parameterized here 
by £ that describes the width of STZ's transition thresh- 
olds distribution. Xoo is the long time limiting value of 
the effective temperature @. Note that the assumption 
that A is small results in a separation of time scales be- 
tween the fast variables A and to and the slow variable x- 
An important feature of the constitutive equations is that 
the onset of homogenous unbounded plastic flow results 
from an exchange of dynamic stability of the bistable 
field to, cf. Eq. (fTl) . For s < 1, i.e. a deviatoric stress 
below s y , the stable fixed-point is to = ±sgn(s), corre- 
sponding to jamming, D pl = 0. On the other hand, for 
s> 1, the stable fixed-point is m=l/s, corresponding to 



flow DP 1 ^ (with A = 



-i/x 



). This observation justifies 



our previous choice of normalizing stresses by the dy- 
namic yield stress s y . The consequences of the exchange 
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FIG. 1: (Color online) The time-dependent radius of the cav- 
ity for asymptotically stable (solid line) and asymptotically 
unstable growth (dashed line). Here and below the parame- 
ters of the model are x°° — 0.13, x(0) = 0.1, fi = 50, p = 1, 
£ = 15. In these runs K — 100. (Inset) m as a function 
of £ at t — 10 4 ro for the two cases. For the bounded growth 
(solid line), m near the free boundary £ = 1 increases towards 
the jamming fixed-point m = l, while in the unstable growth 
(dashed line) m is at the flowing fixed-point m = 1/s with 
s > 1. Wherever 1 there exists a plastic flow, while where 
m — 1 the material is jammed. 



of stability in the homogeneous equations is further dis- 
cussed below in relation to the cavitation instability in 
our highly inhomogeneous configuration. 

The full set of equations was solved numerically. To 
avoid dealing with an infinite time-dependent domain we 
applied the following time-dependent coordinate trans- 
formation 



€ = R(t)/r 



(15) 



This transformation allows us to integrate the equations 
in the time-independent finite domain £ € [0,1]. The 
price is that new terms are generated in the equations. 
This domain was discretized uniformly with 301 nodes. 
Controlling the equations at small distance required the 
introduction of an artificial viscosity on the RHS of Eq. 
([7]). The term introduced is pifJ^v, with r) chosen of the 
order of the square of space discretization over the time 
discretization. Time and length are measured in units of 
To and R(0) respectively. A and m are set initially (£ = 0) 
to their respective fixed-points. 

The first result to be discussed is the existence of a 
cavitation instability, such that for a°° below a threshold 
value a th the circular hole attains an asymptotic stable 
radius, whereas for a°° > a th the growth of the cavity is 
unbounded, signaling a catastrophic failure. The respec- 
tive dynamics are demonstrated in Fig. [TJ 

The value of o- th can be estimated as shown in the ap- 
pendix of Ref. by considering the quasi-static and 
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incompressible analog of the present problem. For the 
parameters at hand (except for K — oo and £ = 10) this 
procedure yielded the estimate a th « 5 in agreement with 
the present results. For ductile materials such stresses are 
hard to achieve and indeed cavitation instabilities were 
observed in the past in ductile materials reinforced by 
or bonded to brittle materials 0- To understand the 
existence of instability note that in typical situations 
the stress at infinity results first in a rapid elastic re- 
laxation to the equilibrium solution of Eq. ([9]); this is 
how the level of applied stress is transmitted to the free 
boundary. When the deviatoric stress exceeds the mate- 
rial yield stress s„, plastic flow initiates, expanding the 
cavity beyond the radius of elastic equilibrium. During 
this process energy is dissipated and the stress tends to 
relax towards s y . When the expansion is slow, the stress 
indeed relaxes to s y , the material becomes jammed and 
the growth is bounded. However, when a°° is sufficiently 
large the rate of growth is such that plastic relaxation 
fails to reduce the level of stress below s y , resulting in a 
cavitation instability. This discussion is highlighted by 
the inset of Fig. [1] where the profiles of m are presented. 

One learns considerable amount of the physics of the 
problem by analyzing the velocity of the free boundary 
in the unstable phase. This velocity has three typical 
regimes; the first is a transient in which s relaxes from its 
initial value <j°° to a typical value somewhat larger than 
Sy. The second regime is very interesting, characterized 
by the fact that the radius R(t) is the only typical scale 
in the problem. We thus expect the velocity of the free 
boundary to satisfy an equation of the form 



R/R = u{a 00 ) > , 



where lu is positive and time-independent in this regime. 
This exponential growth is accompanied with a self- 
similar solution in all the fields of the problem, which 
can depend only on £, cf. Eq. (|15|) . As an example 
we present in Fig. [2] the deviatoric stress s(£) and the 
effective temperature x(?) f° r different times, including 
the onset of the third regime in which self-similar solu- 
tions break down due to inertial effects. This last regime 
manifests itself when the velocity of the free boundary be- 
comes a finite fraction of the typical elastic wave speed; 
when this happens we have a new length scale that can 
be made from the typical wave speed times time, break- 
ing the self similarity. In order to estimate the velocity in 
which inertia becomes important we note that at low ve- 
locities almost all the elastic energy release goes to plastic 
dissipation (below we show that the external work at in- 
finity is negligible). At high velocities we expect that a 
finite fraction of the elastic energy release translates into 
kinetic energy. We thus want to compare the kinetic en- 
ergy density of the medium with the plastic dissipation 
density. The typical kinetic energy density near the free 
boundary is given by pR 2 /2. The plastic dissipation den- 
sity is given by sD^Toe 1 ^ . Consulting Eq. (TIT))) , one 
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FIG. 2: (Color online) The deviatoric stress as a function of 
the variable £ at different times. Note the transient short 
time, the regime of self similarity and the break down of self- 
similarity at long times. (Inset) The effective temperature 
x(£). Note the indistinguishable curves in the self-similar 
regime. In fact, even after the breakdown of self-similarity 
X does not change significantly. In the plastically deforming 
region s > 1 and \ > x(0) (plastic deformation produces disor- 
der). Here the plastic region extends over R(t) <r <6.25R(t) 
(0.16< 1). 
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FIG. 3: (Color online) The velocity of the cavity as function 
of time, for two different values of the bulk modulus K. 



observes that near the free boundary D^TQe 1 ^ is of the 
order of unity and s>s y . Therefore the ratio of the two 
energy densities is of the order of unity when 



R 



hjh 1 



(17) 



Thus in this problem inertia becomes important at 
smaller velocities than in the context of brittle fracture 
where the estimate is v ~ c s [8] (recall that J s y / fi is 
significantly smaller than unity for most materials). 
In Fig. [3] we present the velocity of the free boundary 
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FIG. 4: (Color online) The rate of change of the various en- 
ergies in the problem as a function of time. 



as a function of time for two values of the bulk mod- 
ulus K. Observe the transient regime, the exponential 
growth regime, where the problem exhibits self-similar 
solutions, and finally the attainment of an asymptotic 
velocity when inertial effects become important. In our 
simulations \J s y j p «0.14, and indeed the cross-over for 
K = 100 from the self-similar to the inertial-dominated 
regime occurs where Eq. (|17p predicts. Increasing K by 
two order of magnitude (i.e. approaching the incompress- 
ible limit) results in a further decrease in the cross-over 
point and in the asymptotic velocity. Qualitatively one 
understands this as the effect of faster longitudinal waves 
whose velocity ~ \jKj p that are able to "steal" energy 
from the process near the free boundary. 

Finally, it is worthwhile to consider the energy ex- 
change mechanisms in this problem. Note that in the 
standard approach of Linear Elasticity Fracture Mechan- 
ics energy balance is used to set up the "equation of mo- 
tion" for a crack [8] . On the contrary here we solved the 
boundary value problem, including the explicit evolution 
of the free boundary, and it is interesting to go back with 
this solution at hand to examine the energy exchange. 
Energy comes in four types: the elastic strain energy 
E e i, the kinetic energy Ek, the plastic dissipation W p 
and the external work W ex t- We are mainly interested 
in the time rate of change of these quantities; these are 
given as the following spatial integrals over the domain 
[R(t),R(t)) (the limit R — > oo is briefly discussed below) 
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Wext = 2wa°°R(t)R(t) , W p i=4% sD pl rdr .(18) 

JR(t) 

Conservation of energy implies that 



The limit R — > oo is not well posed as it results in various 
divergent behaviors. We expect the leading term of v and 
v at oo to be ~ r _1 (even for a finite K). In that case 
both Ek and Ek diverge. Moreover, E e i diverges as well 
for a finite K since p approaches a constant at oo. These 
are all peculiarities of 2D that should not affect the basic 
physics we describe. 

In Fig. [4] we plot the different rates of change of energy 
appearing in Eqs. (fT8|) . integrated up to 300i?(i). W ext 
equals the sum of the three other time derivatives such 
that Eq. (fTi?)) is satisfied. Note also that W ex t is sig- 
nificantly smaller than the other three time derivatives. 
Under this condition, the Fig. [4] demonstrates explicitly 
how the release of elastic energy drives the "crack" , con- 
verting this energy into plastic dissipation and kinetic 
energy. 

One question beyond the scope of this Letter is the 
shape stability of the perfectly symmetric cavity. This 
issue will be picked-up in a forthcoming publication. In 
addition, the STZ approach should be applied to cracks 
in less symmetric situations. This requires significant 
investment in tensorial generalizations and in boundary 
tracking algorithms, which are yet other issues under 
study. 
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